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Introduction 

Over  the  last  decade  the  nonlinear  interaction  of  disturbances  in  supersonic  boundary  lay¬ 
ers  has  been  extensively  investigated.  It  has  become  possible  due  to  accumulated  experience  in 
theoretical  and  experimental  studies  of  nonlinear  disturbance  evolution  in  subsonic  flows. 

In  the  case  of  compressible  gas  flows,  especially  at  supersonic  velocities,  the  theoretical  re¬ 
searches  of  nonlinear  evolution  are  strongly  complicated  by  the  following  circumstances.  Consid¬ 
eration  of  temperature  and  density  perturbations  raises  the  order  of  differential  equations,  which 
results  in  an  increase  of  the  calculation  volume.  At  high  Mach  numbers,  additional  unstable  modes 
appear.  According  to  the  linear  theory,  the  most  unstable  first  mode  waves  are  three-dimensional 
(3D)  or  oblique,  the  direction  of  their  fronts  does  not  coincide  with  the  basic  flow  direction. 

For  explanation  of  dynamics  of  unstable  waves  observed  at  introduction  to  a  supersonic 
boundary  layer  of  inspected  disturbances  of  large  intensity,  the  nonlinear  models  of  interaction 
of  disturbances  are  considered.  In  the  models  of  weaklynonlinear  interactions,  mostly  elaborated 
for  subsonic  flows,  one  distinguishes  a  coupling  mechanism  in  the  resonant  triads  and  the  com¬ 
binational  interactions  of  two  oscillations,  comprising  the  processes  of  waves  self-action  (self¬ 
oscillation  modes).  The  possibilities  of  realization  of  three-waves  interaction  with  the  order 
0(s  2)  for  supersonic  flows  are  discussed  in  [1], 

In  the  present  work  the  study  of  the  action  in  the  resonant  triads  including  the  stationary 
modes  and  also  the  study  of  the  second  non-linear  mechanism  realization  possibility  -  interac¬ 
tion  in  pairs  -  is  conducted.  Historically,  it  was  just  this,  and  more  specifically,  the  study  of  the 
finite  amplitude  waves  self-action  effects  was  the  beginning  of  the  near-wall  nonlinearity  simu¬ 
lation. 

The  method  and  results  of  numerical  simulation  of  development  of  disturbances  in  a  su¬ 
personic  boundary  layer  on  a  flat  plate  at  high  Reynolds  numbers  are  presented  in  the  paper.  The 
solution  is  built  by  expansion  with  respect  to  the  small  parameter.  The  contributions  of  linear 
and  quadratic  terms  are  taken  into  account.  The  resonance  interaction  of  the  fundamental 
(2D)mode  of  linear  instability  with  (3D) mode  and  stationary  (3D) mode  in  the  parametric  region 
is  studied.  The  results  are  compared  with  experimental  data  under  controlled  conditions.  The 
results  obtained  show  that  the  suggested  model  accounts  successfully  for  experimental  results 
and  that  the  resonance  mechanism  can  play  an  important  role  in  formation  of  the  spectrum  of 
dominating  disturbances  in  a  supersonic  boundary  layer. 

In  the  present  paper,  only  the  interaction  of  hydrodynamic  waves  which  exponentially  van¬ 
ishes  at  infinity  is  considered,  i.e.,  the  problems  similar  to  the  case  of  low-speed  flows  are  dis¬ 
cussed. 

Formulation  of  the  problem 

The  governing  equations  for  disturbance  evolution  in  a  supersonic  boundary  layer  are  the 
well  known  Navier  -  Stokes  equations  of  continuity,  energy  and  state. 

The  flow  parameters  can  be  presented  as  a  sum 

Q(x,  y,  z,  t)=Qb(x,  y,  z)+  £  Z(x,  y,  z,  t) 
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where  Qb(x,  y,  z)  is  the  solution  to  the  stationary  equations  of  motion,  s  Z  is  the  disturbances 

(s  »!)• 

We  consider  the  evolution  of  disturbances  in  supersonic  boundary  layer  over  flat  plate  at 
high  Reynolds  numbers  Re  =  (Upc/v  e)1/2,  where  Ue,  ve  —  free-stream  velocity  and  kinematic 
viscosity,  x  -  distance  from  the  leading  edge  to  a  reference  point. 

Therefore  parallel  flows  is  a  good  aproximation  for  the  basic  flow. 

The  disturbances  can  be  described  by  a  system  of  nonlinear  equations  which  depend  on  the 
basic  flow  parameters.  For  weak  nonlinearity  it  is  possible  to  restrict  the  solution  by  the  contri¬ 
bution  from  linear  and  quadratic  terms.  Besides  we  shall  take  into  account  viscosity  and  termal 
conductivity  only  in  the  linear  terms  in  highest  derivative.  Then  the  eight-component  vector- 
function  Z (u,  uy,  v,  p  ,T,  Tr,  w,  wy)  can  be  presented  in  the  operational  notation  as  a  solution  of 
the  system  of  differential  equations 

LZ=£-  Mlq^qki).  (1) 

Here  u,  v,  w  are  velocity  perturbations  in  x,  y,  z-directions,  T,  p  are  perturbations  of  temperature 
and  pressure,  q  is  a  four-component  vector,  subscripts  y,  t,  x,  z  mean  the  transport  coefficients 
(viscosity  and  heat  conductivity). 

The  solution  of  (1)  can  be  obtained  by  using  an  expansion  on  the  small  parameter  £  and  a 
two-scale  expansion  on  the  x  -coordinate,  i.e.,  we  introduce  «fast»  (x{)  and  «slow»  scales  X=  £ x. 
A  possibility  of  introducing  two  scales  is  justified  by  a  large  difference  in  phase  and  amplitude 
changing  rates.  Thus, 

d/dx=d/dxi  +s  d/dX,  Z  =  Z°+sZ>. 

In  this  case  Z  satisfied  following  equation  and  boundary  conditions: 

L0  Z °=0,  z/=  z3°=z5°=  $7  °=0  aty  =  0,  y  =  oo .  (2) 

The  zeroth-order  solution  Z  can  be  taken  in  the  form: 

Z °=Real(Z  cij(X) Z\X,  y)exp(i  0  ,)).  (3) 

Coefficients  aj,  the  phases  ©  ;=  I  a/x)dx+  ft  f—cOjt  and  ZI,J  depend  on  slow  variable  X.  Solu¬ 
tion  Z0j  is  governed  by  equation  (A,B,C,D,G)  are  matrixes  which  depend  on  flow  parameters): 
f-i  co  jA +/  a  jB+i  P  jC+G)  Zoj  +D d  Z°>/dy=0 

with  boundary  conditions  (2).  Equation  (3)  subject  to  homogeneous  boundary  conditions  (2) 
form  eigenvalue  problem  which  provides  complex  value  a  ,  for  given  ( co p ,),  Real(  a ,)  being 
the  streamwise  wave  number  and  -/?«( ap  being  the  growth  rate  of  primary  wave.  In  the  next 
order  we  obtain  the  following  set  of  coupled  equations  (mark  «asterisk»  means  complex  conju¬ 
gate): 

L°  Z1  =  -  Z  BZoj  exp(;  ®})).(da/dX)  +M (z0mk,  z0n%), 

where  vector  M  represents  nonlinear  terms.  Applying  usual  solvability  conditions  (orthogonal¬ 
ity  to  the  solution  W0j  of  adjoint  problem)  yields  the  set  of  amplitude  equations : 

da/dX=  -Im(  cq  )cij  +kIa2a3QX p(  ), 

da2/dX~  -Im( a.2  z)a2  +k2a 2a3  exp(  ;A^2 )» 
da3/dX=  —Im{  a 3  )a3  +k3a2a2  exp(  iA<p 3 ), 
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kj=  1  (Mj  W 0])dy/\  (B z0J '  W °>)dy, 

where  A <p  -  phases  which  are  small  for  resonant  triplet. 

We  are  interested,  first  of  all,  in  evolution  of  the  (2D)  linear  unstable  mode  with  (3D)mode 
and  stationary  (3D)mode  which  satisfy  the  resonant  conditions: 

P  i=  P  2+  P  3>  ®  co  2+  co  3. 

For  frame  of  the  combinational  interactions  we  come  to  the  following  set,  describing  the 
spatial  behavior  of  wave  amplitudes: 

da1/dX=  \-Im(  a  /)+(£;,/  _/ 1 ai\+Ei  2,-2  \a2\pP1  'Vii, 

dcil/dX—  [— 7i lll(  (X  ^(E  2,2,-2  2\~  ^ P 2. 1  .-1  \ctp')P  1  ] Cl2> 

Pj=\  W°j dL^/dcCjdy,  j=l,2. 

Here  E  are  the  non-linear  coefficients,  connecting  the  secondary  harmonics  and  original  waves 
(3).  Let  us  write  down  the  structure  of  several  of  them: 

EIX.,=  \  i  _i,Z°i)  W *dy  +  \  MjCZ^  U,Z°E  W 0>  dy, 

Eu,-2=\  Mj(Zvt2_2,Z°i)  W °'dy+\ M^Z*  ly2,Z°.2)  W 0]dy+j  M^Z*  i,_2,Z°2)  W 0]dy. 

Here  Zvl  are  the  secondary  waves,  emerging  by  non-linear  self-  and  cross-  interactions;  Zvt1_1  and 
Zvt2_  2  are  the  zero  harmonics,  describing  deformations  of  the  mean  flow  characteristics;  Z2\  j  is 
the  overtone,  being  created  by  the  first  wave,  its  phase  is  ©  =2  ©  ,;  ZvlL2  is  the  summary  secon¬ 
dary  wave  with  the  phase  0  =  ©  y+  0  2;  Z'Y  2  is  the  differential  secondary  wave  with  the  phase 
©=©;-©,. 


Numerical  results  and  discussion 

In  this  work  the  total  temperature  was  kept  constant  and  equal  to  310  K,  the  Prandtl  num¬ 
ber  was  Pr  =  0.72,  specific  heat  ratio  was  y  =1.4,  the  parameters  corresponding  to  the  wind 
experiments  of  [2]. 

The  nonlinear  interaction  for  resonant  wave  triads  in  the  supersonic  boundary  layer  on  a 
flat  plate  were  numerically  investigated  but  we  will  not  represent  these. 

In  the  present  paper  the  situation  arising  at  introduction  to  a  boundary  layer  on  a  slice  has 
exposed  to  the  theoretical  analysis  at  a  Mach  number  Me=2  of  inspected  disturbances  the  large 
intensity  suffices  [2].  The  authors  have  called  observable  dynamics  of  disturbances  as  "abnor¬ 
mal".  As  against  earlier  considered,  she  is  characterized  by  a  number  of  features  to  explain 
which  within  the  framework  of  used  non-linear  models  it  fails. 

At  "abnormal"  dynamics  is  established,  that  2D  disturbances  are  most  unstable.  In  an  ini¬ 
tial  spectrum  two  wave  packets  on  multiple  frequencies  (f  =  1 0  kHz  -  subharmonic  frequency 
are  arrested  and  f2  =  20  kHz  -  base  frequency),  the  packet  dominates  on  f,  the  maximum  of 
intensity  is  necessary  on  2D  wave.  Downstream  mainly  the  preference  (2D)  character  of  wave 
spectrums  does  not  change,  the  strong  growth  of  disturbances  on  frequency  f  is  watched  (almost 
in  10  times),  the  amplification  of  disturbances  on  a  base  frequency/)  is  a  little  bit  less,  but  all 
the  same  it  is  much  more  linear.  It  results  that  the  laminar-turbulent  transition  happens  on  20% 
closer  to  a  leading  edge.  There  is  the  distortion  of  a  mean  stationary  velocity  U  crosswise  of 
boundary  layer,  which  width  was  increased.  The  phase  velocities  of  disturbances  were  above 
linear  on  30-40%. 
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Fig.  1  The  linear  increments  of  waves  (1  —  in  fj.  2  —  in  fj). 


Fig.  2.  The  secondary  zero  harmonics  m,-i  (1) 
and  U2 -2  (2),  mean  profile  U  (dash-line). 


It  was  clear  that  the  observable  features  are  called  by  non-linear  character  of  development 
of  disturbances  distinct  from  three-wave  resonant.  In  the  present  paper  only  the  combinational 
interactions  of  two  (2D)  oscillations  of  is  considered. 

In  [2]  the  source  of  inspected  disturbances  placed  at  Re  =  497,  and  the  measurements  are 
conducted  in  range  624  <  Re  <  846.  For  these  parameters  calculations  also  are  conducted.  In 
Fig.  1  we  are  exhibited  increments  at  of  linear  waves  on  a  subharmonic  frequency  /)  and  base 
frequency  f2  (digit  1  and  2).  The  initial  position  of  a  subharmonic  at  Re  =  624  -  near  to  the 
lower  branch  of  a  neutral  curve,  the  linear  increment  is  augmented  with  growth  Re,  not  reaching 
a  maximum  in  final  cross-section  of  measurements.  The  basic  wave  places  near  to  a  maximum 
of  a  linear  increment,  down-stream  its  increment  decreases. 

In  Fig.  2  we  present  the  secondary  harmonicses  u2  ,  and  u2_2  at  initial  Re  on  a  background 
of  a  mean  stationary  profile  U  (shaped  line).  Full  deformation  U  will  be  recorded  A  U=\a2\2  Uj 
i+\a2\2  u2,-2  also  will  be  determined  by  amplitudes  of  incident  waves.  The  registration  of 
nonlinearity  is  visible  that  results  in  the  greater  property  of  being  filled  up  of  a  profile  in  wall- 
rigion  and  occurrence  of  a  defect  at  an  outer  boundary,  that  augments  a  boundary  layer  thick¬ 
ness.  It  completely  correlates  with  experimental  observations. 

The  direction  of  the  influence  of  non-linear  processes  of  amplitudes  primary  waves  it  is 
possible  to  analyze  considering  non-linear  factors.  Let's  unite  factors  which  are  in  charge  of  self¬ 
effect  and  for  a  cross  interaction.  In  Fig.  3  they  are  exhibited  for  the  researched  range  Re.  The 
positive  values  of  the  non-linear  members 
will  call  accessory  (concerning  linear) 
growth  of  amplitudes  and  negative  will 
put  to  its  decrease.  From  a  figure  it  be¬ 
comes  clear  that  self-effect  of  a  subhar¬ 
monic  £)j  (full  curve  1)  will  put  to 

increase  of  amplitude  aj  and  on  the  con¬ 
trary  -  self-effect  of  a  wave  on  f2  will 
moderate  its  amplitude,  this  process  will 
rise  with  increase  Re  (shaped  curve  1). 

The  influence  of  a  base  wave  on  a 
subharmonic  (Ei  2  _2)  consists  of  increas¬ 
ing  amplitude  by  last  at  small  Re  and  of 
decreasing  downstream  (full  curve  2).  For 
a  base  wave  (f2)  availability  in  a  spectrum 


Fig.  3.  Non-linear  coefficients  for  the  amplitude  equations. 
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of  a  subharmonic  always  results  in  growth 
a2  (shaped  curve  2). 

It  is  necessary  to  illustrate  consid¬ 
ered  model  on  behaviour  of  amplitudes 
(Fig.  4).  In  calculation  initial  values  a2 
and  a?  is  corresponded  to  experimental 
total  intensities  of  initial  wave  packets, 
and  7/72=5,  and  ai/a2=2. 

In  figure  shaped  lines  exhibit  linear 
dynamics  of  amplitudes,  appropriate  to 
their  linear  increments,  solid  -  are  shown 
non-linear  amplitudes  at  combinative 
process,  and  by  dash-and-dot  lines  are 
given  non-linear  amplitudes  in  a  mode  of 
self-effect. 

For  a  subharmonic  the  nonlinearity 
always  results  in  heightened  growth  of 
amplitude.  Thus  a  linear  degree  of  growth 
(relation  of  final  and  initial  amplitudes)  is  determined  by  value  ~  3  and  value  -  ~  6  for  non¬ 
linear  thus  at  the  expense  of  nonlinearity  the  amplitude  can  bring  up  almost  in  2  times.  The  basic 
contribution  introduces  to  process  self-effect,  registration  of  presence  of  a  wave  on  basic  to  fre¬ 
quency  lowers  a  level  of  accessory  amplification  because  of  transfer  to  her  of  a  part  of  energy  a 
little. 

For  this  base  wave  because  of  considerable  reduction  linear  increments  on  a  researched  in¬ 
terval  Re  the  fading  linear  is  watched  the  amplitudes  a2,  nonlinearity  in  a  mode  of  self-effect 
practically  has  not  an  effect  on  its  value  and  the  influence  of  a  subharmonic  results  to  small  to 
growth  of  this  amplitude.  A  low-level  level  of  influence  of  nonlinearity  is  determined  by  a 
smallness  of  an  initial  amplitude  of  this  wave. 

In  non-linear  process  there  is  also  phase  progression  of  phases  which  can  result  in  to  change 
of  wave  numbers  Real(a)  and  phase  velocities  of  disturbances  c=  co  /Real(  a ).  In  [2]  is  marked 
that  experimentally  particular  values  of  phase  velocities  c  ~  0.7-0.72  on  the  average  on  30- 
40%exceeded  appropriate  values  of  own  waves  supersonic  boundary  layer  (c~  0.52-0.55).  The 
phase  velocities  of  waves  at  non-linear  stage  are  determined.  The  values  of  phase  velocities  with 
non-linear  by  the  corrections  become  very  close  to  experimentally  particular  (c  ~  0.69-0.75). 

Summarizing  said  we  shall  mark  that  within  the  framework  of  the  considered  combinative 
model  qualitatively  prominent  feature  of  dynamics  of  inspected  disturbances  of  heightened  in¬ 
tensity  -  considerable  signal  amplification  as  contrasted  to  linear,  the  distortion  of  mean  velocity 
diagrams  correctly  are  described  in  the  field  of  outer  boundaries  resulting  in  to  growth  the 
thickness  of  the  boundary  layer  and  increase  of  phase  velocities  of  waves. 

This  work  was  supported  by  the  Russian  Foundation  for  Basic  Research  (grant  No.  00-01- 
00828). 
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